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Resumen 

Larios-Tlali, H v Torres-Benites, E., Quevedo-Nolascol, A., 
Martmez-Menes, M. R., & Salgado-Transito, J. A. (mayo- 
junio, 2015). Riesgo de inundacion en la subcuenca del no 
La Antigua, Veracruz, Mexico. Tecnologi'a y Ciencias del Agua, 
6(3), 39-56. 

El objetivo del presente trabajo fue aplicar el modelo de 
simulacion hidraulica HEC-RAS para delimitar las zonas 
con riesgo a inundacion en la subcuenca del rio La Antigua, 
Veracruz, Mexico, utilizando la metodologfa del Centro 
Nacional de Prevencion de Desastres (Cenapred), que 
considera la vulnerabilidad y el peligro de un sistema ante la 
ocurrencia de un fenomeno extremo perturbador. El modelo 
HEC-RAS usa ecuaciones diferenciales determimsticas que 
simulan el comportamiento y la dinamica de los niveles 
de agua en secciones transversales del cauce; su aplicacion 
bajo condiciones de flujo no permanente y regimen mixto 
permitio definir las areas de inundacion alcanzadas 
durante una avenida extraordinaria. Para la simulacion 
hidraulica se construyeron hidrogramas unitarios sinteticos 
con datos de precipitacion registrados durante un evento 
hidrometeorologico extremo; la geometrfa del cauce se 
obtuvo con la extension HEC-GeoRAS, herramienta que 
extrajo y proceso informacion georreferenciada del Modelo 
Digital de Elevacion (MDE) LiDAR tipo terreno y los 
coeficientes de rugosidad (q de Manning) se estimaron de 
acuerdo con las caractensticas ffsicas del canal. Para calibrar 
los datos estimados por el modelo se utilizo informacion 
hidrometrica (niveles registrados) de la estacion 28003 
Cardel. La aplicacion de tres estadfsticas de prueba: R, RMSE 
y DMA demostraron un alto ajuste entre datos observados 
versus simulados, respaldando la aplicabilidad del modelo 
HEC-RAS para simular, predecir y delimitar zonas de riesgo 
a inundacion. 

Palabras clave: simulacion hidraulica, inundacion, 

hidrograma, fenomenos hidrometeorologicos extremos, 
calibracion, riesgo a inundacion. 


Abstract 

Larios-Tlali, H, Torres-Benites, E., Quevedo-Nolascol, A, 
Martmez-Menes, M. R., & Salgado-Transito, J. A. (May-June, 
2015). Flood Risk in a Sub-Basin of the La Antigua River, Veracruz. 
Mexico. Water Technology and Sciences (in Spanish), 6(3), 39- 
56. 

The aim of this study was to apply the HEC-RAS hydraulic 
simulation model to delineate flood hazard areas in a sub-basin of the 
La Antigua River, Veracruz, using the methodology by the National 
Center for Disaster Prevention (Cenapred). This methodology takes 
into account the vulnerability and risk of a system to an extreme 
disturbance. The HEC-RAS model uses deterministic differential 
equations to simulate the behavior and dynamics of water levels 
in cross-sections of channels. By applying it under unsteady and 
mixed flow conditions, areas of flooding resulting from an unusually 
large flood can be defined. For hydraulic simulation, synthetic 
unit hydrographs were constructed using data recorded during 
an extreme hydrometeorological precipitation event. Channel 
geometry was obtained with HEC-GeoRAS, a tool which extracts 
and processes georeferenced information from a LiDAR land 
surface Digital Elevation Model (DEM). Roughness coefficients (q 
Manning) were estimated according to the physical characteristics 
of the channel. To calibrate the estimates calculated by the model, 
hydrometric information (recorded levels) was used from the Cardel 
28003 station. The application of three test statistics— R, RMSE 
and DMA— showed a good fit between simulated and observed 
data, indicating the applicability of the HEC-RAS model to simulate, 
predict and delineate flood hazard zones. 

Keywords: Hydraulic simulation, flood hydrograph, extreme 
hydrometeorological events, calibration, flood. 
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Introduction 

Mexico, por su ubicacion geografica, es un pais 
vulnerable a la ocurrencia de diversos fenome- 
nos hidrometeorologicos, principalmente a la 
formacion y desarrollo de ciclones tropicales 
(Jauregui, 2003) en la costa del Pacffico y del 
Golfo de Mexico. Un efecto de estos fenomenos 
-en cualquiera de sus etapas- son las inunda- 
ciones originadas por lluvias torrenciales; sin 
embargo, no es el unico factor que determina 
la magnitud de estas catastrofes, tambien 
incluye la distribucion espacial de la lluvia, la 
morfometria de los cauces, la infraestructura, la 
topografia, el relieve, la cobertura vegetal, el uso 
de suelo y las unidades edafologicas (Cenapred, 
2001). 

La zona costera del estado de Veracruz, 
por su ubicacion geografica, orograffa, clima 
y asentamientos humanos en zonas de alto 
riesgo, tambien es un area vulnerable a estos 
fenomenos. En las dos ultimas decadas (de 1990 
a 2010), los municipios de Tuxpan, Tamiahua, 
Tecolutla, Nautla, San Andres Tuxtla y La An- 
tigua sufrieron los efectos de cinco tormentas 
tropicales y diez huracanes (GEV, 2011). 

El paso cercano de tormentas o ciclones 
tropicales pueden ser causas de inundaciones, 
pero tambien lo pueden ser las depresiones 
tropicales menores e incluso nortes (Acevedo & 
Luna-Dfaz, 2006). Entre 1950 y 1994 se registra- 
ron 475 inundaciones en la region Golfo Centro, 
de las cuales 255 ocurrieron en la region Panuco 
(CSVA, 2007), la mayoria de estas catastrofes 
se relacionan con la ocurrencia de fenomenos 
hidrometeorologicos extremos. 

El potencial hfdrico superficial del estado 
de Veracruz es de los mas altos de Mexico, su 
escurrimiento superficial anual medio se esti- 
ma en 121 millones de metros cubicos, lo que 
representa 33% del escurrimiento superficial de 
todo el pafs (CSVA, 2005). Veracruz tiene una su- 
perficie de 72 815 km 2 , 73% de su territorio esta 
por debajo de los 200 msnm (Tejeda, 2006). El 
gobierno del estado de Veracruz (GEV, 2011) es- 
tima una superficie inundable de 6 275 km 2 (8% 
del territorio estatal), donde habitan 1 252 251 


personas en 2 189 localidades de 118 municipios, 
lo que indica que uno de cada seis veracruzanos 
se ve afectado por las inundaciones. 

Por lo anterior, es necesario e imprescindible 
realizar un estudio integral para estas zonas con 
alto riesgo a inundacion, el cual debe proponer 
acciones, obras y programas que prevengan, 
reduzcan, atiendan y mitiguen los danos ocasio- 
nados por tales catastrofes. Una de las medidas 
no estructurales para disminuir el riesgo por 
inundacion es la delimitacion de dichas zonas. 

El presente trabajo ha combinado el uso 
de sistemas de informacion geografica (SIG), 
tecnicas de percepcion remota, trabajos de 
prospeccion en campo, analisis de informacion 
hidrometeorologica y aplicacion de un modelo 
matematico de simulacion hidraulica ( HEC - 
RAS , version 4.1.0) para delimitar las zonas 
con riesgo a inundacion en la subcuenca del rio 
La Antigua, Veracruz, una de las regiones del 
estado mas danadas por inundaciones en las 
ultimos anos. 

Los resultados obtenidos en la simulacion 
hidraulica fueron comparados con los niveles 
registrados en la estacion hidrometrica 28003 
Cardel (durante la ocurrencia de un fenomeno 
hidrometeorologico extremo), para ello se em- 
plearon tres estadfsticas de prueba (R, RMSE y 
DMA), cuyos valores obtenidos demostraron 
un alto ajuste entre datos observados versus 
simulados. 

Area de estudio 

La cuenca del rio La Antigua se ubica entre los 
paralelos 19° 05' y 19° 34' latitud norte y los 
meridianos -96° 06' y -97° 16' longitud oeste. 
Tiene un area aproximada de 2 827 km 2 , que se 
distribuye en los estados de Puebla y Veracruz. 
El rio La Antigua nace en la Sierra Madre Orien- 
tal, con el nombre de rio Sumidero, al oriente 
de la localidad Gonzalez Ortega, en el estado 
de Puebla, a una altitud de 3 345 msnm (Espi- 
noza, 2009). La cuenca pertenece a la Region 
Hidrologica 28 Papaloapan (RH28) y a la Su- 
bregion Jamapa y Otros, y esta conformada por 
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cuatro subcuencas: RH28Bc, RH28Bd, RH28Be 
y RH28Bf (figura 1). 

El area de estudio incluye solo la subcuenca 
RH28Be (subcuenca del rio La Antigua), ya que 
en esta zona se registran los mayores danos 
provocados por el desbordamiento del no La 
Antigua. 

La subcuenca RH28Be tiene una superficie 
de 508.3 km 2 , se enmarca dentro de los paralelos 
19° 12' y 19° 23' y los meridianos -96° 21 ' y -97° 
00'. Se distribuye horizontalmente con direccion 
oeste-este en el estado de Veracruz, sobre las 
subprovincias fisiograficas denominadas "Chi- 
conquico" y "Llanura Costera Veracruzana", las 
cuales pertenecen, respectivamente, a la Provin- 
cia Fisiografica X Eje Neovolcanico Transversal 
y a la Provincia Fisiografica XIII Llanura Costera 
del Golfo Sur. 

El clima predominante, segun la clasificacion 
de Koppen, modificada por Enriqueta Garcia, 


es el grupo A (calidos humedos tropicales), 
donde domina la unidad climatica "Awo", con 
una temperatura que oscila entre 18.7° y 25.2 °C; 
en los meses de abril a mayo se presentan los 
valores mas altos y de noviembre a febrero las 
temperaturas mas frias; la temperatura media 
anual es de 22.3 °C. Por otro lado, la precipita- 
tion media anual es de 1 393 mm; en la parte 
alta de la subcuenca (zona montanosa) se han 
registrado valores que superan los 2 000 mm, y 
en la parte central se registran valores entre 900 
y 1 000 mm. La evaporation media anual es de 
1 340 mm; en los meses de marzo a mayo es mas 
intensa y en el mes de diciembre se presentan 
los valores mfnimos (SMN, 2008). De acuerdo 
con el conjunto de datos vectoriales de uso de 
suelo y vegetacion serie IV, escala 1:250 000 del 
INEGI (2012), en la subcuenca RH28Be existen 
seis usos de suelo y vegetacion: agricultura, 
pastizal cultivado, bosque de encino, bosque de 
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Figura 1. Ubicacion geografica de la cuenca del rio La Antigua, con base en la Red Hidrografica escala 1:50 000 del Instituto 

Nacional de Estadfstica y Geograffa (INEGI), 2010. 
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comferas, selva baja caducifolia y selva mediana 
subcaducifolia; el uso de suelo predominante 
es la agricultura, ocupa una superficie de 322.3 
km 2 , lo que representa 63.4% del area total. Por 
otro lado, con base en el Sistema de Clasificacion 
de Suelos WRB-2000, la unidad edafologica 
de mayor distribution es el vertisol pelico, al 
ocupar una superficie de 271.7 km 2 , lo que re- 
presenta 53%, seguido del acrisol ortico y litosol, 
con 110.7 y 98.3 km 2 , respectivamente. 

Materiales y metodos 

Description del ntodelo 

El HEC-RAS (Hydrologic Engineering Center- 
River Analisis System) es un modelo nume- 
rico creado por el Cuerpo de Ingenieros de la 
Armada de los Estados Unidos, de uso libre 
y en continuo desarrollo (Turegano, 2012). El 
procedimiento basico de calculo de este software 
se fundamenta en la solution de la ecuacion de 
energia, donde las perdidas de energia por fric- 
tion se estiman usando la ecuacion de Manning; 
las perdidas locales por contraction y expansion 
del flujo se calculan como una fraction del cam- 
bio en la velocidad entre dos secciones (USACE, 
2002). Este programa utiliza el metodo de paso 
estandar (Standard Step Method), el cual consis- 
te en determinar el tirante en una seccion a partir 
del tirante asociado con el escurrimiento de una 
seccion adyacente (aguas arriba o aguas abajo, 
dependiendo del regimen de escurrimiento). El 
balance de energia entre estas dos secciones se 
hace resolviendo la siguiente ecuacion: 

y 2 +z 2+ ^=r 1+ z l+ ^ + h e a) 
2 g 2 g 

donde Y y Y 2 = tirante en la seccion trasversal; 
Z v Z 2 = elevaciones del fondo del canal prin- 
cipal; V v V 2 = velocidades promedio; a v a 2 = 
coeficientes de velocidad; g = aceleracion de la 
gravedad, y h g = perdida de energia (Ceballos, 
2011). 

Segun Timbadiya, Patel y Porey (2011), este 
razonamiento simplifica el flujo de agua como 


un flujo unidimensional, asume que el tirante y 
velocidad solo varian en direction longitudinal 
del canal, cuyo eje se supone aproximadamente 
a una linea recta con velocidad constante en 
cualquier punto de la seccion trasversal. 

Delimitation de las unidades de analisis 

La delimitation de la unidad de analisis esta- 
blece las interacciones entre los componentes 
del sistema y los procesos que en el ocurren 
durante la simulation hidraulica. El mapa base 
se elaboro utilizando los datos vectoriales de la 
RH28B, extraida de la Red Hidrologica escala 
1:50 000 reportada por INEGI en su pagina web 
(consultado el 1° de julio de 2013, disponible en 
http:/ / www.inegi.org.mx). 

El analisis, como se ha senalado, incluyo solo 
la subcuenca RH28Be, debido a que en esta zona 
se han registrado los may ores danos provocados 
por el desbordamiento del rio La Antigua; sin 
embargo, considerando que el escurrimiento 
registrado por la estacion hidrometrica 28003 es 
el acumulado de RH28Be y RH28Bf (figura 2), 
cuyos datos fueron utilizados en la calibration, 
se modelo el comportamiento y la dinamica de 
la red hidrologica RH28Be, tomando en cuenta 
el escurrimiento generado por RH28Bf. 

En la figura 3 se observa la suma del escu- 
rrimiento de RH28Bf al tramo "Sta. Maria", a 
traves del afluente conocido como "Pescados". 
De esta forma se modelaron ocho cauces de 
RH28Be y uno de RH28Bf (ver figura 2). 

Para cada afluente se requirio de informa- 
cion especifica de su zona de dominio, por lo 
que se dividio al area de estudio en nueve uni- 
dades de respuesta hidrologica (URH); cada una 
corresponde a un afluente o tramo considerado 
en la modelacion. Los parametros de las URH 
que se utilizaron para estimar los caudales se 
muestran en el cuadro 1. 

Estimation de los hidrogramas de entrada 

En la subcuenca RH28Be, los registros de pre- 
cipitacion son mas abundantes que los de escu- 
rrimiento, por ello se eligio un metodo indirecto 
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Figura 2. Distribucion de las subcuencas y unidades de respuesta hidrologica (URH). Para delimitar las URH se empleo WMS, 

ver. 8.3 y el MDE LiDAR tipo terreno del INEGI, 2013. 



Figura 3. Confluencia del rio Pescados (subcuenca RH28Bf) y el rio La Antigua (subcuenca RH28Be). 


para construir un hidrograma de escorrentia 
superficial, asociado con una lluvia de duracion 
y profundidad unitaria. El metodo utilizado 


fue el hidrograma adimensional curvilmeo del 
SCS (Soil Conservation Service), que puede ser 
representado por un hidrograma unitario trian- 
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Cuadro 1. Parametros morfometricos de las URH. 


URH 

Parametro 

Nombre 

Num. 

Longitud del 
cauce (km) 

Cota maxima 
(msnm) 

Cota minima 
(msnm) 

Superficie 

(km 2 ) 

Pendiente (%) 

Xochilapa 

I 

38.8 

1556 

260 

193.96 

0.032 

Actipan 

II 

16.3 

984 

260 

44.91 

0.044 

La Bomba 

III 

35.8 

900 

182 

53.96 

0.020 

Huilapa 

IV 

28.3 

722 

182 

40.57 

0.019 

Pescados 

V 

113.9 

3 721 

98 

1 631.67 

0.031 

Sta. Maria 

VI 

63.0 

1556 

98 

314.02 

0.023 

La Bomba 2 

VII 

46.4 

900 

73 

110.06 

0.018 

Sta. Maria 2 

VIII 

120.0 

3 721 

73 

1 950.06 

0.030 

Parte Baja 

IX 

142.2 

3 721 

5 

2 145.15 

0.026 


gular (HUT). El metodo propuesto por Mockus 
(1957) genera un gasto pico en funcion del area 
de la cuenca y el tiempo pico segun la siguiente 
ecuacion: 

0.555hp e A (2) 

Del analisis de varios hidrogramas, Mockus 
concluye que el tiempo base y el tiempo pico se 
relacionan mediante la expresion: 


h = 2 - 67 t v 

(3) 

el tiempo pico se expresa: 


d 

K= — + t r 

v 2 r 

(4) 

VO 

o 

ii 

s-J*" 

(5) 

El tiempo de concentracion ( t ) 
estimar con la ecuacion de Kirpich: 

se puede 

t c = 0.000325 TL 

(6) 

La duracion del exceso de lluvia (d) se 
calcula para cuencas grandes con la siguiente 
ecuacion: 

d e =2jT c 

(7) 


donde q^= gasto pico unitario (m 3 /s/mm); 


hp e = altura de precipitacion efectiva (mm); 
A = area de la cuenca (km 2 ); t p = tiempo pico 
(h); t b = tiempo base (h); t r = tiempo retraso (h); 
L = longitud del cauce principal (m), y S = pen- 
diente del cauce (%). 

Una vez generado el HUT con la metodolo- 
gia de Mockus, se estimo el hidrograma adimen- 
sional curvilmeo del SCS, que suaviza la forma 
del HUT; el hidrograma del SCS conserva las 
mismas unidades de tiempo y volumen que el 
HUT. Para obtener el hidrograma unitario del 
SCS se multiplican las ordenadas por el q v y las 
abscisas por el t (Barrios & Olaya, 2007). 

El hidrograma unitario se construyo con los 
datos de precipitacion registrados en 25 esta- 
ciones meteorologicas (ubicadas dentro y en la 
periferia de la subcuenca RH28Be) durante la 
ocurrencia de un fenomeno hidrometeorologico 
extremo que impacto las costas de Veracruz el 17 
de septiembre 2010: el huracan Karl. Los datos 
de precipitacion fueron proporcionados por la 
Direction Tecnica del Organismo de Cuenca 
Golfo Centro (OCGC) de la Comision Nacional 
del Agua (Conagua). Con la information de las 
estaciones se calculo la precipitacion media pon- 
derada de cada URH utilizando el metodo de las 
isoyetas, el cual establece un area de influencia 
delimitada por lmeas que unen los puntos de 
igual valor de precipitacion de acuerdo con la 
siguiente expresion: 
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donde D = altura de precipitacion media (mm); 
a. = area entre cada dos isoyetas (km 2 ); D = pro- 
medio de precipitacion entre dos isoyetas (mm), 
y A = area total de la cuenca (km 2 ). 

En la figura 4 se muestran los hidrogramas 
adimensionales curvilmeos del SCS estimados 
para las nueve URH, y en el cuadro 2 se tienen 
los valores de t y Q, valores que construyen 
dichos hidrogramas. 

Geometria de cauces 

La geometria de los cauces se obtuvo modelan- 
do con HEC-GeoRAS , una extension del HEC- 
RAS que extrae la informacion georreferenciada 
de un MDE. En este trabajo se utilizo el MDE 
LiDAR tipo terreno, con una resolution de 5 m 
en la horizontal y 0.25 m en la vertical. 

Se initio digitalizando el centro del cauce 
principal de la red hidrologica RH28Be, comen- 
zando de aguas arriba hacia aguas abajo; los 
centros de cada afluente tributario se unieron 


conforme avanzo la digitalization. Posterior- 
mente se dibujaron las margenes, primero 
la margen izquierda (del cauce principal) y 
despues su margen derecha; de igual forma se 
unieron las margenes de los afluentes tributarios 
de aguas arriba hacia agua abajo. La longitud de 
los centros de cauce digitalizados suman 186.5 
km y las margenes suman 374.4 km (figura 5). 

Para generar las secciones transversales se 
tomaron en cuenta las siguientes condiciones: a) 
el trazo de la seccion parte de la margen izquier- 
da hacia la margen derecha; b) las secciones se 
trazan de aguas arriba hacia aguas abajo y son 
perpendiculares al centro del cauce; c) las sec- 
ciones no se pueden interceptar entre ellas, y d) 
las secciones transversales solo pueden cruzar 
una sola vez al centro y margenes del cauce. En 
total se crearon 2 435 secciones transversales a lo 
largo de los nueve cauces o tramos que integran 
la red hidrologica digitalizada, el intervalo entre 
cada seccion fue de 20 m y la longitud de trazo 
(o ancho de seccion) oscilo entre 100 y 1 200 m. 
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Figura 4. Hidrogramas adimensionales curvilmeos del SCS considerados en la simulacion hidraulica. 
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Secciones 

transversales 



Figura 5. Secciones transversales de los rios La Bomba, Huilapa y La Bomba 2, obtenidas con la extension HEC-GeoRAS. 


Determination del coeficiente de Manning 

El coeficiente de rugosidad de Manning repre- 
senta la resistencia al flujo de agua en cauces 
sobre las llanuras de inundacion; el valor es 
mas alto cuando existe mayor rugosidad en la 
superficie de contacto con el flujo (Kumar, 2013). 
Este coeficiente deberfa ser calibrado a partir 
de information obtenida en campo (perfiles 
y secciones de nos); sin embargo, cuando no 
se dispone de datos de aforo, se recomienda 
utilizar gufas — para la selection del coeficien- 
te — que propongan valores para cauces de 
caracterfsticas similares o valores obtenidos de 
datos experimentales (Turegano, 2012). 

El modelo requiere la asignacion de un coefi- 
ciente de rugosidad para la margen izquierda, 
derecha y centro de cada seccion transversal. 
El valor se establecio utilizando la "Tabla del 
coeficiente de rugosidad de Manning", pro- 
puesta por Chow, Maidment y Mays (1994), la 
cual establece el coeficiente en relation con las 


caracterfsticas y condiciones mas sobresalientes 
del canal, como rugosidad del fondo del cauce, 
cobertura y tipo de vegetation, y sinuosidad 
del canal, entre otras. La asignacion de este 
coeficiente se apoyo en puntos de verification 
en campo efectuados sobre el no La Antigua y 
sus afluentes. 

En la figura 6a se observa que el fondo del 
cauce posee distintas caracterfsticas de rugo- 
sidad; por lo tanto, el fondo de cada seccion 
transversal posee un coeficiente de Manning 
acorde con las condiciones del sitio. Tambien se 
observa que ambas margenes pueden ser muy 
diferentes en una misma seccion (figura 6b), por 
ello cada margen posee un coeficiente propio. 
En total se obtuvieron 7 305 valores, estos datos 
se ingresaron en una base de datos dentro del 
entorno de HEC-RAS. 

Ademas de los coeficientes de Manning 
correspondientes a cada seccion transversal, 
el modelo requiere: a) la separation entre el 
extremo izquierdo de una seccion y el extremo 
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a) Asignacion de coeficiente de Manning para el fondo del cauce 





Planicie de inundacion 
sin vegetacion 


Planicie de inundacion con 
vegetacion exuberante 


b) Asignacion de coeficiente de Manning para las margenes y llanuras de inundacion 


Figura 6. Asignacion de los valores de Manning segun tablas de Chow, Maidment y Mays (1994). 


izquierdo de la seccion de abajo; b) la distancia 
que existe entre el cruce de una seccion al eje y 
el cruce de la seccion de abajo, y c) la separacion 
entre el extremo derecho de una seccion y el 
extremo derecho de la seccion de abajo. Estas 
distancias fueron medidas y se ingresaron en 
una base de datos dentro del entorno de HEC- 
RAS. 

Condiciones defrontera y regimen deflujo 

Una vez capturada la base de datos, se ingresa- 
ron los hidrogramas unitarios del SCS estima- 
dos para cada URH. Dada la sinuosidad de los 
afluentes se eligio simular bajo un regimen de 
flujo mixto, el cual consider a las condiciones de 
frontera aguas arriba y abajo del cauce. 


En la confluencia de un afluente tributario 
y el cauce principal no se requirio establecer 
condiciones entre dichas uniones, ya que el 
modelo considera automaticamente condiciones 
de contorno interior, solo se definieron las con- 
diciones en los extremos aguas arriba y abajo; 
para ello se ingresaron las pendientes a la que 
se asocio un tirante uniforme al inicio y al final 
de cada afluente. 

Indices de evaluacion estadistica de ajuste del 
modelo 

La calibration busca la correlation entre los 
valores de las variables generadas por el mo- 
delo y las medidas en el sistema real (Ponce, 
1986). Moriasi (2007) recomienda el uso de los 
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siguientes indices estadisticos de prueba para 
determinar la correlacion entre ambos valores: 
La raiz del error cuadratico medio (RMSE) 
que se calcula con la siguiente ecuacion: 


RMSE = 


2U*,df} 


N 


(9) 


La desviacion media acumulada (DMA), 
segun Molnar (2011), esta dada por la siguiente 
ecuacion: 

= Ll — il (io) 

N 

Y el coeficiente de correlacion (. R ) dado por: 


R = - 






( 11 ) 


donde a. es el dato estimado por el modelo; t., el 
dato que se asume como medido u observado; 
N, el mimero de observaciones o estimaciones; a, 
el promedio de los datos estimados por el mode- 
lo; y t, el promedio de los datos que se asumen 
como medidos (Cai, Liu, Lei, & Pereira, 2007; 
Cervantes, Arteaga, Vazquez, Ojeda, & Queve- 
do, 2013). En los dos primeros indices, el valor 0 
corresponde a un ajuste perfecto, mientras que 
un valor mas alto representa un menor ajuste. 
Para R, el valor 1.0 indica una alta correlacion; 
entre 0.5 y < 1.0, una buena correlacion; entre 0.0 
y < 0.5, una baja correlacion, y 0 que no existe 
correlacion (ASCE, 1993). 

Resultados 

A1 usar el modelo HEC-RAS se simulo el com- 
portamiento hidraulico de los afluentes que 
integran la red hidrologica de la subcuenca 
RH28Be, esto como respuesta a un gasto o escu- 
rrimiento que generaron las lluvias del huracan 
Karl el 17 de septiembre de 2010. 

La simulacion hidraulica permitio obtener el 
tirante o nivel del flujo de agua en cada una de 


las secciones transversales, sin embargo, estos 
resultados no pueden ser considerados como 
"aceptables" si no se realiza una prueba de ca- 
libracion. Segun Pappenberger, Beven, Horritt y 
Blazkova (2005), la importancia de la calibracion 
radica en la posibilidad de aplicar sus predic- 
ciones para describir, explicar y comprender un 
sistema real. 

Para llevar a cabo este proceso, se efectuo 
antes un analisis de sensibilidad que consistio 
en realizar diferentes simulaciones variando 
dos parametros del modelo: a) distanciamiento 
entre secciones transversales y b) coeficiente de 
rugosidad. 

Los resultados obtenidos por las diferentes 
simulaciones (en especffico niveles o tirantes) 
se compararon con los niveles observados en la 
estacion hidrometrica 28003 (entre las 00:00 y las 
23:00 horas del dfa 18 de septiembre de 2010); 
para ello se utilizaron tres indices estadisticos 
de prueba que determinaron la correlacion entre 
ambos valores. 

La corrida o simulacion hidraulica elegida 
para los fines del presente estudio fue aquella 
en la que se obtuvo el mejor ajuste — al aplicar 
los tres indices de evaluacion estadfstica — entre 
valores observados versus simulados (cuadro 3). 

Con la informacion del cuadro 3 se ajusto un 
modelo lineal entre los niveles medidos y obser- 
vados que se pesentan en la figura 7, donde hay 
una menor relacion cuando el nivel supera los 
11 m; el coeficiente de determinacion (R 2 ) para 
el conjunto de datos es de 0.896. 

Las tecnicas graficas de evaluacion son herra- 
mientas que facilitan la comparacion visual de 
los datos simulados versus observados y permi- 
ten una valoracion inmediata de los resultados 
(ASCE, 1993); sin embargo, los metodos graficos 
revelan poca informacion sobre la calidad de los 
resultados y el funcionamiento del modelo. Para 
evaluar la capacidad de HEC-RAS se emplearon 
tres indices estadisticos de prueba: RMSE, DMA 
y R. El cuadro 4 muestra el resultado de estos 
indices utilizados en la prueba de calibracion. 

Los valores obtenidos indican una buena 
correspondencia entre datos simulados y 
observados, esto confirma que la modelacion 
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Cuadro 3. Niveles observados versus simulados (18 de septiembre de 2010). 
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Figura 7. Niveles observados en la estacion hidrometrica 28003 Cardel y simulados por HEC-RAS. 


Cuadro 4. Resultados del analisis de comparacion de niveles observados y simulados por el modelo HEC-RAS. 


Indice de evaluacion estadistica 

Valor 

Raiz del error cuadratico medio 

0.866 

Desviacion media acumulada 

0.744 

Coeficiente de correlacion 

0.946 
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hidraulica con HEC-RAS es aplicable para de- 
limitar las zonas con riesgo a inundacion de la 
subcuenca RH28Be. 

Delimitation de las zonas con riesgo a 
inundation 

Quaas (2004) menciona que el riesgo depende de 
dos factores fundamentales: el peligro y la vul- 
nerabilidad. Guevara, Quaas y Fernandez (2004) 
definen al peligro como la probabilidad de que 
ocurra un fenomeno potencialmente danino (o 
factor perturbador) durante un cierto periodo 
de tiempo y en un sitio dado. Por otra parte, 
Cardona (2001) define vulnerabilidad como 
la susceptibilidad o propension de un sistema 
expuesto a un fenomeno o factor perturbador. 
El Cenapred (2004) define al factor perturbador 
como un fenomeno natural o generado por el 
hombre que ocasiona danos a un sistema en un 
grado tal que constituye un desastre. 


Establecimiento del factor perturbador 

Para determinar el peligro se establecio como 
factor perturbador el area indudable (o llanura 
de inundacion) obtenida a traves de la simu- 
lacion con HEC-RAS. En la figura 8 se observa 
que el factor perturbador no representa un 
peligro para todas las comunidades, esto se 
debe principalmente al relieve de la subcuen- 
ca y la distancia que existe entre el centro de 
poblacion y el rio La Antigua o afluentes. La 
probabilidad de que ocurra un fenomeno o 
evento potencialmente danino para una po- 
blacion no depende solo del comportamiento 
hidraulico del cauce. 

De las 96 comunidades que existen en la 
subcuenca RH28Be (donde habitan 34 877 per- 
sonas, segun el Censo de Poblacion y Vivienda 
2010), solo para 23 de ellas representa un peli- 
gro el desbordamiento del rio La Antigua (ver 
cuadro 5), fenomeno representado por el nivel 



Figura 8. Factor perturbador y peligro en un tramo del rio La Antigua. 
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o tirante que el agua alcanzo en las secciones 
trasversales durante la simulation hidraulica. 

La frecuencia o recurrencia de inundaciones 
depende de la frecuencia de precipitaciones 
excepcionalmente fuertes; un estudio realiza- 
do por la Conagua (2011) indica que el gasto 
generado por las lluvias del huracan Karl tiene 
un periodo de retorno de 500 anos, por lo cual, 
segun Ahmed (2008), su frecuencia situa a la 
mayorfa de las comunidades (de RH28Be) en un 
nivel bajo de amenaza. Esta forma de evaluation 
se usa por lo comun en pafses como Australia, 
Estados Unidos y Suiza; en el presente trabajo se 
utilizo la metodologfa del Cenapred (2011), que 
aplica como criterios fundamentales el peligro 
y la vulnerabilidad. 

Determination del indice de vulnerabilidad 

Segun Cardona (2001), la vulnerabilidad ff- 
sica corresponde al analisis de las viviendas 
para conocer que tan susceptibles son a una 
inundacion. El analisis se hizo con base en la 
tipologfa de las viviendas en el ambito local; 
una vivienda es mas vulnerable si los mate- 
riales de construccion son de baja resistencia, 
como adobe, carton o lamina; y sera menos 
vulnerable si esta hecha de cemento, ladrillo 
o tabique (materiales de mayor resistencia). 
La information referente al tipo de material 
utilizado para su construccion se extrajo del 
Censo de Poblacion y Vivienda 2010, editado 
por el INEGI (2011). 

El fndice de vulnerabilidad se determino so- 
lo para las 23 comunidades a las que representa 
un peligro el desbordamiento del rfo La Antigua 
y sus afluentes (como resultado de la simula- 
cion hidraulica). De esta forma se obtuvo que 
11 comunidades, donde habitan 988 personas, 
poseen un fndice de vulnerabilidad alto y muy 
alto; cinco comunidades, donde habitan 4 737 
personas, tienen un fndice de vulnerabilidad 
medio; y siete comunidades, donde habitan 28 
043 personas, poseen un fndice de vulnerabili- 
dad bajo y muy bajo. 


Cartografia de riesgo a inundation 

La base conceptual para elaborar el mapa de 
riesgo a inundacion se describe con mas detalle 
en Cenapred (2011); sin embargo, se exponen a 
continuation los criterios mas sobresalientes de 
esta metodologfa, que derivo en la cartografia 
de riesgo: a) riesgo alto es toda zona cuyo tirante 
es > 1 m y la velocidad es > 1 m/s; b) riesgo 
medio es la zona cuyo tirante esta comprendido 
entre 0.4 y 1 m, y la velocidad fluctua entre 0.4 
y 1 m/s; c) riesgo bajo es la zona cuyo tirante es 
< 0.4 m y la velocidad es < 0.4 m/ s. 

Con base en estos criterios y los resultados 
obtenidos en la modelacion hidraulica, se esta- 
blecio el nivel de riesgo de las 23 comunidades 
(cuadro 5) y se genero un mapa de riesgo a 
inundacion. La figura 9 representa un bosquejo 
de este mapa correspondiente a la Parte Baja del 
rio La Antigua (URH IX). El mapa se edito a una 
escala tal que permite observar la distribution 
de las comunidades, el peligro y la vulnerabili- 
dad de estas a una inundacion. 

La modelacion permitio demostrar que las 
lluvias asociadas con el huracan Karl provo- 
caron el desbordamiento de los afluentes que 
integran la red hidrologica de RH28Be; sin 
embargo, no todos los sitios afectados repre- 
sentan un riesgo directo para los habitantes de 
las 23 comunidades, debido a que cinco de los 
nueve tramos simulados impactan solo en zonas 
agrfcolas, pastizales, o relictos de bosque y selva 
(los tramos contenidos en las URH II, III, IV, V y 
VII). Este hecho no debe exceptuar la atencion 
de dichas zonas, pues segun los reportes de 
Valdez y Arreola (2011), en Veracruz — a causa 
del huracan Karl — se perdio casi 100% de las 
cosechas en cultivos como la papaya y el mafz, 
lo que en terminos monetarios represento mas 
de 100 millones de pesos (mdp), y en el sector 
pecuario se perdieron mas de 1 200 cabezas de 
ganado bovino, 11 000 cerdos y 20 000 aves. 

Por otro lado, las URH I, VI, VIII y IX ad- 
quieren mayor importancia social. En ellas se 
localiza 100% de las localidades con algun nivel 
de riesgo a inundacion. Los reportes de Armen- 
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Cuadro 5. Comunidades de la subcuenca RH28Be bajo algun nivel de riesgo. 


Num. 

Comunidad 

Municipio 

Poblacion total 

Latitud 

Longitud 

Nivel de riesgo 

1 

Jose Cardel 

La Antigua 

24 846 

19.37 

-96.38 

Alto 

2 

Juan Yobal Zerna 

Puente Nacional 

4 

19.29 

-96.61 

Alto 

3 

Los Cuates 

Puente Nacional 

4 

19.28 

-96.63 

Alto 

4 

Puente Nacional 

Puente Nacional 

347 

19.32 

-96.48 

Alto 

5 

Salmoral 

La Antigua 

753 

19.34 

-96.35 

Alto 

6 

Salmoral parte del Sur 

Paso de Ovejas 

1 161 

19.34 

-96.36 

Alto 

7 

San Jose de Enmedio 

Puente Nacional 

6 

19.30 

-96.58 

Alto 

8 

San Jose de La Azud 

Puente Nacional 

9 

19.30 

-96.58 

Alto 

9 

Yucatan 

Paso de Ovejas 

543 

19.34 

-96.36 

Alto 

10 

Colalco 

Puente Nacional 

9 

19.33 

96.50 

Medio 

11 

El Coyolar 

Puente Nacional 

390 

19.29 

-96.61 

Medio 

12 

El Cuajilote 

Puente Nacional 

258 

19.30 

-96.58 

Medio 

13 

Mata Grande 

Paso de Ovejas 

345 

19.35 

-96.39 

Medio 

14 

Paso de la Cruz 

Puente Nacional 

7 

19.35 

-96.39 

Medio 

15 

San Jose Chipila 

Puente Nacional 

402 

19.29 

-96.59 

Medio 

16 

Santa Maria Tatetla 

Jalcomulco 

1777 

19.27 

-96.71 

Medio 

17 

Santo Nino de Atocha 

Puente Nacional 

5 

19.30 

-96.66 

Medio 

18 

Chichicaxtle 

Puente Nacional 

1910 

19.34 

-96.47 

Bajo 

19 

El Cedro 

Puente Nacional 

71 

19.30 

-96.57 

Bajo 

20 

El Manguito 

Paso de Ovejas 

173 

19.36 

-96.38 

Bajo 

21 

Paso Marino 

Puente Nacional 

388 

19.33 

-96.46 

Bajo 

22 

Tomas Vallejo 

Puente Nacional 

3 

19.29 

96.59 

Bajo 

23 

Xopilapa 

Tenampa 

357 

19.28 

-96.78 

Bajo 


96° 28' 96° 26' 96° 24' 96° 22' N 



Figura 9. Mapas de riesgo a inundacion correspondiente al tramo Parte Baja (URH IX) de la RH28Be. 
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ta et al. (2011) mencionan que el huracan dejo 
un saldo aproximado de 500 000 damnificados 
y danos estimado en 70 mil mdp. 

Conclusiones 

Los datos de precipitacion registrados durante 
el huracan Karl permitieron estimar los hidro- 
gramas de escorrentia superficial asociados 
con dicho evento meteorologico. Asimismo, 
los puntos de verificacion en campo de uso de 
suelo y vegetacion facilitaron la asignacion del 
valor T| de Manning (o coeficiente de rugosidad) 
mediante el uso de tablas con valores preesta- 
blecidos. El MDE de alta resolution permitio 
representar las condiciones morfometricas 
del no La Antigua y afluentes secundarios 
sin requerir un levantamiento topografico, 
aunque esto ultimo no significa que el trabajo 
de topografia sea innecesario o deba olvidarse. 
Toda esta information procesada para su uso en 
HEC-RAS permitio delimitar zonas de riesgo a 
inundacion. 

Los resultados obtenidos por el modelo se 
compararon con datos observados en la estacion 
hidrometrica 28003 mediante el uso de indices 
estadisticos de prueba, los cuales indicaron una 
alta correlation entre ambas series de datos, 
respaldando de esta forma la aplicabilidad del 
modelo HEC-RAS para predecir la magnitud de 
una inundacion, catastrofe muy frecuente en la 
subcuenca RH28Be. 

Segun los resultados obtenidos en la mo- 
delacion hidraulica del no La Antigua y sus 
afluentes, asi como el uso de la metodologia 
del Cenapred (2011) para elaborar mapas de 
riesgo a inundacion (que aplica como criterios 
fundamentales el peligro y la vulnerabilidad), 
se obtuvo que nueve localidades, donde habitan 
27 673 personas, poseen un nivel de riesgo alto; 
ocho un nivel medio con 3 193 personas, y seis 
en un nivel bajo con 2 902 personas. Estos datos 
indican que 96.8% de la poblacion total que 
habita en la subcuenca RH28Be radica en una 
zona con riesgo a inundacion. 

De la elaboration y analisis del mapa de 
riesgo a inundacion se concluye que las comu- 


nidades mas cercanas al rio La Antigua y cuyas 
viviendas en su mayoria son de materiales de 
baja resistencia, son las de mayor riesgo. Por el 
contrario, las localidades mas alejadas al rio La 
Antigua y cuyas viviendas son de materiales 
resistentes, se clasifican en localidades de bajo 
riesgo. 

Cualquiera que sea el caso, desbordamientos 
en zonas agricolas o centros de poblacion, se re- 
comienda la implementation de medidas estruc- 
turales y no estructurales contra inundaciones, a 
fin de prevenir, controlar y/o mitigar los danos 
ocasionados por estos eventos meteorologicos 
extremos. 

Referencias 

Acevedo, R. F., & Luna-Dfaz, R A. (2006). Principales 
fenomenos meteorologicos que afectaron al estado de Veracruz 
en el ano 2005. Inundaciones en el Estado de Veracruz 2005 
(pp. 39-49). Veracruz, Mexico: Universidad Veracruzana. 
Ahmed, H. K. M. (2008). Aplicaciones del modelo HEC-RAS 
para el analisis delflujo no permanente con superficie libre (290 
pp.). Quito, Ecuador: Escuela Politecnica Nacional. 

ASCE (1993). Criteria for evaluation of watershed models. 
American Society of Civil Engineers. Journal Irrigation 
Drainage Eng., 119(3), 429-442. 

Armenta, M. S., Lopez, J. C., Rodriguez, E., Allan, E., Del 
Amo, S., Gomez, A., Mac Swiney, M. C., Niembro, A., 
Sanchez, O., Vazquez, M., & Velazquez, N. (2011). La 
restauracion ecologica como estrategia para la reduction del 
riesgo de desastres ante inundaciones. Inundaciones 2010 en 
el estado de Veracruz (pp. 196- 236). Mexico, DF: CVICDT. 
Barrios, P. M. I., & Olaya, M. E. J. (2007). Calculo y analisis 
de hidrogramas para el flujo torrencial del 22 de junio 
de 2006 ocurrido en la microcuenca "El Salto", Ibague, 
Colombia. Avances en Recursos Hidrdulicos, 16, 31-40. 

Cai, J., Liu, Y., Lei, T., & Pereira, L. S. (2007). Estimating 
reference evapotranspiration with the FAO Penman- 
Monteith equation using daily weather forecast messages. 
Agricultural and Forest Meteorology, 145, 22-35. 

Cardona, A. O. D. (2001). La necesidad de repensar de manera 
hoh'stica los conceptos de vulnerabilidad y riesgo. Una cn'tica 
y una revision necesaria para la gestion (18 pp.). Ciudad 
de Panama: Red de Estudios Sociales en Prevencion de 
Desastres en America Latina. 

Ceballos, L. J. D. (2011). Modelacion hidraulica y 
morfodinamica de cauces sinuosos aplicacion a la 
Quebrada La Marinilla (ANT). Boletm de Ciencias de la 
Tierra, 30, 107-118. 

Cenapred (2001). Diagnostico de peligros e identification de 
riesgos de desastres en Mexico (232 pp.). Serie Atlas Nacional 


K VlSojOUDdJ^ 


Larios-Tlali dal., Riesgo de inundacion en la subcuenca del rio La Antigua, Veracruz, Mexico 


de Riesgos de la Republica Mexicana. Centro Nacional de 
Prevencion de Desastres. 

Cenapred (2004). Gui'a Bdsica para la elaboradon de Atlas 
Estatales y Municipales de Peligros y Riesgo (389 pp.). Serie 
Atlas Nacional de Riesgos de la Republica Mexicana. 
Centro Nacional de Prevencion de Desastres. 

Cenapred (2011). Metodologi'a para la elaboradon de mapas de 
Riesgo por inundaciones en zonas urbanas (73 pp.). Serie 
Atlas Nacional de Riesgos de la Republica Mexicana. 
Centro Nacional de Prevencion de Desastres. 

Cervantes, O. R., Arteaga, R., Vazquez, M. A., Ojeda, W., 
&. Quevedo, A. (2013). Comparacion de modelos para 
estimar la presion real de vapor de agua. Tecnologi'a y 
Ciendas del Agua , 4(2), 37-54. 

Chow, V. T., Maidment, D. R., & Mays, W. (1994). Hidrologi'a 
aplicada (584 pp.). Santafe de Bogota: McGraw Hill. 

Conagua (2011). Estudio Integral del no La Antigua del estado 
de Veracruz (pp. 145-146). Mexico, DF: Comision Nacional 
del Agua. 

CSVA (2005). Programa hidrdulico estatal. Resumen ejecutivo (40 
pp.) Veracruz, Mexico: Consejo del Sistema Veracruzano 
del Agua. 

CSVA (2007). Balance Hidrdulico de la Cuenca del Rio Sedeno (27 
pp.). Veracruz, Mexico: Consejo del Sistema Veracruzano 
del Agua. 

Espinoza, T. B. (2009). Evaluadon de los recursos hi'dricos de 
la cuenca del no La Antigua, Veracruz: estado actual y ante 
escenarios de cambio climdtico (102 pp.). Veracruz, Mexico: 
Universidad Veracruzana. 

GEV (2011). Plan Veracruzano de Desarrollo 2011-1016 (303 
pp.). Veracruz, Mexico: Gobierno del Estado de Veracruz. 

Guevara, E., Quaas, R., & Fernandez, G. (2004). Gui'a bdsica 
para la elaboradon de atlas estatales y municipales de peligros 
y riesgos: lineamientos generates para la elaboradon de atlas 
de riesgos (pp. 15-29). Mexico, DF: Sistema Nacional de 
Proteccion Civil-Cenapred. 

INEGI (2011). Censo de Poblacion y Vivienda 2010. Mexico, DF: 
Instituto Nacional de Estadistica y Geografia. 

INEGI (2012). Conjunto de Datos Vectoriales de Uso de Suelo y 
Vegetacion serie IV, escala 1:250 000 (132 pp.). Mexico, DF: 
Instituto Nacional de Estadistica y Geografia. 

Jauregui, E. (2003). Climatology of Land Falling Hurricanes 
and Tropical Storms in Mexico. Atmosfera, 16(A), 193-204. 

Kumar, P. P. (2013). HEC-RAS Model for Mannnig-s 
Roughness: A Case Study. Scientific Research, 3(3), 97-101. 

Mockus, V. (1957). Use of Storm and Watershed Characteristics in 
Synthetic Hydrograph Analysis and Application. Sacramento, 
USA: American Geophysical Union, Pacific Southwest 
Region. 

Molnar, P. (2011). Calibration. Watershed Modelling, SS 2011. 
Zurich: Institute of Environmental Engineering, Chair of 
Hydrology and Water Resources Management, ETH. 


Moriasi, D. N. (2007). Model Evaluation Guidelines for 
Systematic Quantification of Accuracy in Watershed 
Simulation. Transactions of the AS ABE, 50(3), 885-900. 

Pappenberger, F., Beven, K., Horritt, M., & Blazkova, 
S. (2005). Uncertainty in the Calibration of Effective 
Roughness Parameters in HEC-RAS using Inundation 
and Downstream Level observations. Journal of Hydrology, 
302, 46-69. 

Ponce, V. M. (1986). Diffusion Wave Modeling of Catchment 
Dynamics. Journal of Hydrologic Engineering, 109( 4), 637- 
638. 

Quaas, R. (2004). Gui'a bdsica para la elaboradon de atlas estatales 
y municipales de peligros y riesgos (pp. 9-10). Mexico, DF: 
Sistema Nacional de Proteccion Civil-Cenapred. 

SMN (2008). Normales climatologicas [en linea]. Consultado 
el 5 de mayo de 2014. Servicio Meteorologico Nacional. 
Recuperado de http://smn.cna.gob.mx/productos/ 
normales/ cnormales.html. 

Tejeda, M. A. (2006). Panoramica de las inundaciones en el 
estado de Veracruz durante 2005. Inundaciones en el estado 
de Veracruz 2005 (pp. 9-20). Veracruz: Universidad 
Veracruzana. 

Timbadiya, V., Patel, P. L., & Porey, P. D. (2011). Calibration 
of HEC-RAS Model on Prediction of Flood for Lower Tapi 
River, India. Journal of Water Resources and Protection, 3(11), 
805-811. 

Turegano, C. F. (2012). Gestion de riesgo de inundacion (254 
pp.). Madrid: Universidad Politecnica de Madrid. 

USACE (2002). HEC-RAS User Manual (25 pp.). Washington, 
DC: US Army Corps of Engineers. 

Valdez, R. V., & Arreola, V. A. (2011). Estrategias para reducir la 
vulnerabilidad de las inundaciones de produccion de tilapia ante 
inundaciones. Inundaciones 2010 en el estado de Veracruz (pp. 
527-548). Veracruz: CVICDT. 

Direccion institucional de los autores 


Ing. Hugo Larios-Tlali 
Dr. Abel Quevedo-Nolasco 
Dr. Mario R. Martinez-Menes 

Colegio de Postgraduados 

Campus Montecillo 

km 36.5 carretera Mexico-Texcoco 

56230 Montecillo, Texcoco, Estado de Mexico, Mexico 

Telefono: +52 (045) 5527 68 4527, (595) 9520 200, 

extensiones 1384 y 1213 

irnr_larios@hotmail . com 

anolasco@colpos.mx 

mmario@colpos .mx 


Tecnologta y 


Ciencias del Agua, vol. VI, num. 3, mayo-junio de 2015 


Ciencias delAgua, vol. VI, num. 3, mayo-junio de 2015 


Larios-Tlali et al, Riesgo de inundation en la subcuenca del rio La Antigua, Veracruz , Mexico 


M. C. Elibeth Torres-Benites 

Manejo Integral de Cuencas, S. A. de C. V. 

Jose Ma. Pino Suarez, Mz. 15, Lt. 19, Col. Netzahualcoyotl 
56234 Texcoco, Estado de Mexico, Mexico 
Telefono: +52 (595) 9529 824 
elyt@micas-cuencas 


Dr. Jorge A. Salgado-Transito 

Colegio Mexicano de Especialistas en Recursos 
Naturales, A. C. 

De las Flores 8, San Luis Huexotla 
Texcoco, Estado de Mexico, Mexico 
Telefono: +52 (595) 9285 236 
arturo .tr ansito@gmail. com 


i( VlSoiOUDdJ^ 



